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Abstract 



a, 

We analyze the combined effect of a Laplacian field and quenched disorder 
for the generation of fractal structures with a study, both numerical and 
theoretical, of the quenched dielectric breakdown model (QDBM). The growth 
dynamics is shown to evolve from the avalanches of invasion percolation (IP) 
to the smooth growth of Laplacian fractals, i. e. diffusion limited aggregation 
— | (DLA) and the dielectric breakdown model (DBM). The fractal dimension 

is strongly reduced with respect to both DBM and IP, due to the combined 

effect of memory and field screening. This implies a specific relation between 
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the fractal dimension of the breakdown structures (dielectric or mechanical) 
and the microscopic properties of disordered materials. 



H " PACS: 61.43.-j; 61.43.Hv; 02.50.+S. 

C3 ' 



The growth of fractal structures is usually described in terms of physical models charac- 
terized by an irreversible dynamics and a degree of self-organization. These models can be 
divided in two broad classes. In the first one the dynamics is stochastic ( modulated by a 
field), and its prominent examples are diffusion limited aggregation (DLA) [[L]], the dielectric 
breakdown model (DBM) and their variants. The second class is characterized by a 
deterministic dynamics in a quenched random medium. An example is Invasion Percolation 
(IP) || and a related model is the Bak and Sneppen (BS) @] model of biological evolution. 
Concerning the growth dynamics, some of these models grow in a smooth way (i.e. DLA 
and DBM), while other ones grow by avalanches (i.e. IP) Q. 

Our understanding of the origin of fractal structures for the two classes is based on the 
idea that an effective screening is present in the scale invariant dynamics Q. However, 
the origin and the properties of this screening are very different. For Laplacian fractals it 
arises from the geometrical screening of the electric field around the structure. For the case 
of quenched dynamics, instead, there is no field and an effective screening develops as a 
memory effect from the quenched dynamics itself. 

The two kinds of dynamical models, stochastic with a field and quenched without a 
field, represent however two extreme limiting cases of real natural phenomena which usu- 
ally present both features. For example dielectric breakdown and fracture propagation in 
disordered solids represent important cases in which these elements are combined |||7|. 

In this letter we address the question of understanding how the combination of these two 
basic elements operates and leads to new phenomena. From a theoretical point of view this 
requires the unification of the concepts developed for the two limiting cases. Together with 
suitable simulations, this allows us to link the microscopic properties of disordered materials 
to the resulting fractal structures. 

In the original DBM the local field E t across a perimeter bond i of the structure is 
related to the probability of growth by pi oc E^ , where 77 is a parameter characterizing the 
strength of the link between the field and the growth probability. The growth probability 
of a bond is then normalized with that of all the other perimeter bonds. If we consider 



the generalization of such a process in a random medium |||7]] each bond can be associated 
to a quenched variable x« G [0,1], extracted from a given probability distribution, which 
defines the strength of the bond (resistivity) with respect to the breakdown. The growth 
process proceeds by breaking the weakest bond. Without the field this means to break the 
bond with mim-fa;,}; this is Invasion Percolation ||. If we introduce the local field Ei, the 
tendency to break will be modulated by the field itself. This requires the introduction of a 
new variable: 

* w = pSf (1) 

which takes into account this effect. Growth at time t occurs at the bond with the smallest 

Vi(t). 

The probability density of x, is usually assumed of a power law form || 

Po (x) = ax 11 - 1 (2) 

where the parameter a (a G [0, oo]) models the strength of the material. Small values of a 
(a << 1) correspond to a fragile material in which most bonds are easy to break. Large 
values of a (a >> 1) correspond instead to a strong material. The limit a — ► oo (77 — > oo) 
implies that growth occurs deterministically at the bond with the largest field. The case 
rj = eliminates the effect of the field and leads to the IP model. On the other hand, there 
is no obvious limit that brings us back to the original (stochastic) DBM. In this respect the 
present model, that we may call Quenched DBM (QDBM), contains the elements of both 
the DBM and IP models, but it has a simple limiting case only towards IP. 

The effect of the field on the quenched variables is to modify the extension of its distri- 
bution for a given bond, as shown in fig. [I]. So, while the extension of x% is on the [0, 1] 
interval for all bonds, the distribution of yi will have an extension depending on the local 
field Ei (fig. [I]). It is convenient to use a normalized field (at time t) 
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where E max (t) is the maximum value of the field among the perimeter bonds at time t. 
Using this field we have yiit) = Xi/lE^t)] 71 , with a probability density (from eq. |2|: 

My) = a[E'M} av y a - 1 (4) 

and therefore the range of the variable yi(t) for the bond with the largest field is always 
< Vi(t) < 1, which is more convenient for the analysis. 

A fundamental characterization of the growth process comes from the acceptation pro- 
files, corresponding to the distributions of the Xi, yi for the grown bonds. In fig. ^| we 
report, for parameter values a = 1 and rj — 1, the distribution of the variables Xi of the 
grown bonds, while the insert refers to the analogous distribution for the y^. 

The behaviour of these distributions as a function of the total growth time t shows the 
self-organization of the system towards specific limiting distributions. The distribution of 
the Xi converges to a function which extends on the entire [0, 1] interval. This implies that a 
value of Xi close to 1 (unfavourable) may actually prevail on all others because of the effect 
of the field. This situation is very different from the case of IP, in which the analogous 
distribution is a theta function with a discontinuity at x = p c < 1, where p c is the critical 
bond percolation probability. The existence of such a threshold implies that the dynamics 
evolves by scale invariant avalanches PUlPl- 

The acceptation profile of the yi (insert of fig. |2|) seems to converge to a theta function, 
but its meaning is quite different than for the case of IP. In fact, the distribution must be 
zero for y > 1 by construction, because at least the variable with the maximum field must 
be smaller than one. Therefore the bonds that are allowed to grow must refer to values of 
y in [0, 1]. On the other hand its approximately flat behaviour is quite non trivial and it 
indicates that the effect of the field dominates with respect to that of the quenched disorder. 
The best candidates for growth are near the tips because of the field, so that memory effects 
are limited with respect to IP and the process is close to deterministic DBM (which would 
give D = 1). 

The fractal structures corresponding to QDBM are shown in fig.||, while in Table | we 



report the values of the QDBM fractal dimensions for different values of a and 77. These 
values give a strong indication for the invariance property 

D f ( m a) = D f (Ti-a) (5) 

This result is quite interesting and non trivial in view of the very different role that the 
parameters r\ and a play in the growth process. It may be explained by observing that if we 
make the transformation 

Xi ^l", ?7 -> // = a/??, (6) 

the ordering of variables yi do not change and the dynamical evolution leads to the same 
fractal cluster. 

From the theoretical side the growth process is deterministic and it arises from the 
quenched disorder modulated by the field of the structure itself. A suitable approach to deal 
with a problem of fractal growth with quenched disorder is the method of the Run Time 
Statistics (RTS) P~|rO|]- It consists in the mapping of a quenched dynamics into a stochastic 



one with cognitive memory. The basic concept of the mapping is the following: if a bond 
has "lost" many times, there is a finite probability that it will never grow, even after an 
infinite time. This introduces an effective dynamical screening which is at the basis of fractal 
properties. The RTS method provides a systematic technique to describe this phenomenon, 
which is typical of models with quenched disorder. 

In the present case the situation is more complex in view of the modulation induced by the 
field. To this purpose it is convenient to introduce the following concepts. To each interface 
bond we assign an effective probability density p T ,t,i{y), where i gives the position of the bond, 
t is the total growth time and r is the time during which the bond has been part of the growth 
interface without being selected ("age"). As soon as a bond becomes part of the growing 
interface (r = 0), its distribution is the original one po,t,i(y) = Po(y) = o\E[{f)\ ar, y a ~ l . After 
a time r the density is modified in a way that is conditional to the growth history of the 



bond. Using the rules of the conditional and composed probabilities |yj and by performing 



an average over quenched disorder, we can derive an expression for the growth probability 
/J> T ,t,i °f the variable i at time t. This implies that i is the extremal variable and the expression 



of its growth probability is 12 



f 1 tt f Ym 

Vr,t,i= dyp Tit4 (y)[[ dy m p TmAm (y m ), (7) 
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where Y m = , E , \ t „ v and the product accounts for the competition between the growing bond 
and the other perimeter bonds. In the same way we can obtain an equation for the update 
from time t to t + 1 of the effective densities: 
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(8) 



gT^r J and we accounted for the change of the electric field of the variable i 
from time t to time t + 1. 

Equations |7|, § provide the mapping of the original quenched problem into a stochastic 
one defined by growth probabilities. The specific realizations of this stochastic process have 
the same statistical weight as those of the original quenched problem. The approximation in 
this mapping consists in neglecting the geometrical correlations between probability densities 
of different bonds, but this can be shown to be exact in the limit of large systems |14] . Other 
approximations of technical nature will be necessary for the practical use of this mapping 
in the FST scheme for the calculation of the fractal dimension. 

From the expression [7| it is possible to derive the invariance property || that we inferred 
from the simulations. In fact one can show that eq. ^ is invariant under the transformation |6| 



|I"2||. We can also study the behaviour of screening effects in QDBM dynamics. An analysis 



of eqs. ^ and || gives |L2| : 

AirA, oc [E' t (t)} av(T+1) (9) 

The electric field of the bond % at time t, expressed in terms of the initial value of the field 



E'iito) (to being the time at which the bond becomes part of the interface) and of the "age" 



r of the bond is, in analogy with DBM ||13|| : 



E[(t = to + t) = E' t (t )e-^T (10) 

where q iT is a function which depends on the growth history. By using eq. 110 into eq. 



one obtains: 

fx T , t>i oc e -«.r«r(^+i) (11) 

This result shows that screening effects in QDBM are stronger than in stochastic DBM 
(where [i T ,t,i oc e~ CT < iT ). This explains therefore why QDBM has a fractal dimension much 
smaller than the usual DBM or DLA. Moreover, the exponential screening of the GPD 
of QDBM, compared to the power law screening (fJ, Tj t,i oc r^rr^ @) corresponding to IP 
and similar models, gives a further evidence for the absence of scale invariant avalanches 
in QDBM. In fact, as discussed also in |I0[ , the shape of the GPD of extremal processes 
determines the shape of the correspondent avalanche size distribution. Therefore, a power 
law distribution for the avalanche sizes can only arise from a power law screening of the 
growth probability of the interface bonds with respect to the growth time. 

Having characterized the quenched dynamics of the QDBM in terms of RTS growth 
probabilities, we can now consider the use of the Fixed Scale Transformation (FST) fl"5H 
method to analyze the fractal properties of the resulting structures. Such a step, however, 
is nontrivial and it requires the following considerations. The original dynamics is quenched 
as in IP, but in the end it does not lead to a critical threshold for the effective probability 
distribution. This implies that we cannot consider the growth by avalanches as in IP, 
which would correspond to a specific implementation of the FST ||10[ . On the contrary, the 
gaussian screening of eq. O implies that a fast convergence of the FST matrix elements will 



be achieved using the same scheme as in DLA and DBM [|T5 |. 

A nontrivial problem, instead, is the identification of the scale invariant dynamics for 
this problem. In IP the simple extremal rules and the avalanche dynamics can be used 



to argue that the small scale dynamics is already scale invariant [|Kj]. In DLA (DBM) the 
situation is more complex and it requires a specific renormalization study which also showed, 
however, that the small scale dynamics is rather close to the scale invariant one |5j. In the 
present case a systematic study of the scale invariant dynamics would require the combined 
renormalization of the effective probability densities and of the electric fields of the bonds. 
At the moment, this appears far too complex. So, inspired by the results of the two cases 
we can treat, we assume that the small scale dynamics is a reasonable approximation to the 
scale invariant one. 

In this way we can proceed to the explicit calculation of the FST matrix elements and 
of the corresponding fractal dimensions. The results are reported in table for different 
values of 77 with a = 1 (the calculation is up to the second order in the FST scheme). 

These results allow us to understand the strong reduction of the fractal dimension for 
QDBM (D = 1.37 for 77 = 1) with respect to both IP (D = 1.8879 |IJ) and DBM (D = 



1.6406 for 1] = 1 and cilinder geometry [15]). On the other hand, the values of D we obtain 
from theory appear to be somewhat larger than those of the simulations. One can speculate 
that this may be due to the approximation used for the scale invariant dynamics. 
We would like to thank G. Caldarelli for useful discussions. 
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FIGURES 
FIG. 1. Schematic picture of the role of a field (Ei) on the dynamics of QDBM. The two bonds 

have the same (flat) distribution for the variable Xi, which represents the microscopic properties 

of the material. However, once the bonds are subjected to different fields (E\ > E2), the effective 

variables for the breakdown process become the y« = Xi/Ei, whose distributions are modulated by 

the field. 

q 

FIG. 2. Time evolution of the acceptation profiles for the variables Xi and j/j (insert) for the 
bonds which grow. For t — > 00 the convergence towards limiting distributions corresponds to the 
self- organized nature of the dynamics. The absence of a threshold value for a(x) implies a smooth 
growth (not by avalanches). 

FIG. 3. Examples of QDBM clusters for different values of a ■ rj. Note the small values of D 
with respect to IP and DBM. 
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TABLES 



7] D(rj;a = l) a D(a;rj = l) 

0.2 1.33 ±0.02 0.2 1.35 ± 0.02 

0.5 1.21 ±0.02 0.5 1.22 ±0.02 

1.0 1.15 ±0.02 1.0 1.15 ±0.02 

2.0 1.07 ±0.01 2.0 1.06 ±0.01 

3.0 1.02 ±0.01 3.0 1.03 ±0.01 
TABLE I. Fractal dimension of QDBM clusters of size 512 x 2048, in cilinder geometry, for 
different values of the parameters a and r/. 

7] D{rj) 

0.2 1.70 

0.5 1.56 

1.0 1.37 

2.0 1.23 

3.0 1.15 
TABLE II. Second order FST computation of the fractal dimension of QDBM, for different 
values of rj with a = 1. 
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